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Starting from the same input as the standard BFKL Pomeron, we directly calculate 
the "hard" Pomeron as a gluonic ladder by using Monte Carlo methods. We reproduce the 
characteristic features of the the BFKL Pomeron and are now also able to evaluate new 
observables. The applicability of the BFKL approach under realistic kinematical conditions 
can be tested and the influence of the running coupling constant examined. 



I. INTRODUCTION 



A large amount of data for reactions at high energies could be explained very succesfully in terms 
of the exchange of the Pomeron. It was introduced into the general framework of reggeon exchanges in an 
ad hoc fashion and provides the dominant contribution. It is therefore a major challenge to understand the 
nature of this Pomeron microscopically in terms of QCD. Most of the data so far involved small momentum 
transfers, and the exchanged Pomeron was referred to as a "soft" Pomeron. In trying to understand this 
phenomenology on the basis of QCD, it can be argued that the underlying mechanism mainly consists of 
the exchange of gluons, is of long range and thus involves non-perturbative dynamics. 

Attempts to understand the Pomeron from a perturbative QCD point of view lead to the concept 
of a "hard" Pomeron that is involved in small distance, high momentum transfer processes. The interest 
in this Pomeron was renewed when recent data from HERA at low x and large Q 2 showed features that 
could be explained with the exchange of such a hard Pomeron, which has a behavior quite different from 
the soft one. The theoretical approach to understand the hard Pomeron goes back to the work of Lipatov 
et al. (BFKL) @,|. These authors were able to infer some properties of the hard Pomeron perturbatively, 
describing it in terms of a reggeized gluon ladder. Their calculation was performed in leading logarithmic 
order (a s logs ~ 1 for a fixed a s < 1 ) where the multiregge kinematics gives the dominant contribution. 

In this paper, we use a numerical method to better understand the features of the hard Pomeron and 
to see to what extent the concept of a perturbative Pomeron as embodied in the BFKL approach remains 
valid. Rather than proceeding as BFKL and solving analytically under the conditions mentioned above, 
where logs goes to infinity, we use a Monte Carlo method to directly evaluate reggeized gluon ladders 
over a broad kinematical range and for different values of a s . This allows us to study the ladder structure 
explicitly and in more detail. In particular, we are interested to see to what extent the kinematics within 
the ladder stays in a region where perturbative QCD can be applied. Our method also allows us to examine 
the consequences of letting the coupling constant run. Our approach differs from previous numerical studies 
which use angular ordering || or start from the BFKL equations [Q. 



II. THE BFKL POMERON 



We will shortly review the approach leading to the BFKL Pomeron. The BFKL Pomeron can be 
viewed as the exchange of gluonic ladders between two particles A and B as shown in Fig.l. The exchanged 
ladder must be a color singlet state. Working in the leading logarithmic approximation, the main contri- 
bution then comes from the multiregge region. In terms of the rapidities yi of the emitted gluons with 
momenta h, in this kinematical region the rapidity intervals, Syi, between two neighboring emitted gluons 
with rapidities y^i and yi, satisfy Syi 3> 1. The BFKL Pomeron thus is a gluon ladder where the rapidities 
along the ladder increase monotonically. The total rapidity interval, Y = Yb — Ya, between the particles A 
and B can be simply expressed as the sum of the individual rapidity intervals: Y — Y^i=i ^Vi ■ 

The main features of the BFKL Pomeron can be discussed by starting from the expression for the 
total inclusive cross section for the reaction A + B — > X, proportional to the imaginary part of the elastic 
amplitude with a gluon ladder exchange; an example with n rungs is shown in Fig.l. The differential cross 
section has the form 

da n (A + B^X) = ^H (y-Y: SyA \T n \ 2 , (1) 

where T n is the amplitude for emission of n hard gluons (gluon jets). The Pomeron exchange consists 
of the sum over all such ladders with different numbers of gluons. The main ingredients in this gluon 
emission amplitude are the non-local, gauge invariant effective gluon vertices, T M , and the reggeized gluon 
propagators. In the cross section, the vertices appear only in a form contracted over color indices and gluon 
polarizations, i.e. as a contraction of two "Lipatov vertices", 

T^(q k ,q k+1 )-T^(q k ,q k+1 ) = A— ±— - — — . (2) 

(q/t± - qfe+i±j + m 

The gluon mass /i is introduced in order to regulate the infrared behaviour. The Lipatov vertices can only be 
used in the case of large rapidity intervals because otherwise higher-order corrections to the vertices become 
important. The reggeization of the gluons is taken into account by replacing the free gluon propagators by 
the reggeized ones: 



e <5j/;e(ti) 



(3) 



ti — /j, 2 U - /j, 2 

where ti = q 2 ~ — q| and e the gluon trajectory. For large momenta e is given by 



The reggeized gluons represent the exchange of color octet gluon ladders. Their use can be interpreted 
as taking into account Sudakov-like suppression; the inclusion of radiative corrections corresponds to the 
probability that there is no emission of additional gluons in rapidity interval Syi. In this way any double 
counting is avoided. Since the infrared divergences coming from the vertices and the regggeized propagators 
cancel each other, one can safely take the limit /i — > as is done in the original BFKL Pomeron. It can be 
shown that the colour octet part of gluon ladder exchanges as in Fig.l indeed yields the reggeized gluon, 
proving the self-consistency of this approach. 

The reaction mechanism that is independent of the properties of the inital particles A and B, the 
exchanged Pomeron, is indicated by the dashed box in Fig.l. For the dashed box one can derive a recursion 
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relation in the number of the emitted gluons which leads to the well-known BFKL equation [|]J|] . For the 
cross section for this subprocess, the interaction of two gluons (with transverse momenta ha± and q_B_i_), 
the BFKL equation predicts in the asymptotic limit 



J-^fyA- ^frfeV^- l , (5) 



where 

A = a s — log2~2.65-a s , B = 14 C(3) ~ 12.6 • a s . (6) 

7T 7T 

It is necessary that the initial and final transverse momenta of the ladder are large such that perturbation 
theory can be expected to hold. This means that the dimensions of the hadrons A and B (Ra and Rb) should 
be small. From Eq.(||) one can now see several characteristic features of the BFKL Pomeron exchange: 



Energy dependence of the total cross section: 

After integration over the gluon transverse momenta, one obtains for a 



' exp(Ay) ~ s A . (7) 



Y 

For a realistic value of the fixed coupling constant, a s = 0.2, one obtains for A = 0.53, which is much 
larger than for the soft Pomeron, while the low-x data indicate A ~ 0.3. So it will be interesting to 
find out what will be the effect of including a running coupling constant in the BFKL Pomeron on 
the intercept. 

• Diffusion pattern of the transverse momentum distribution: 

From Eq. (ra) one also immediately recognizes the diffusion behavior in the logarithm of the transverse 
momenta . For a gluon with transverse momentum q j_ and with rapidity Ya + y the logarithm 
of the transverse momentum will under the condition 1 C i/ C 7 fluctuate around its initial value 
log(q^j_) according to 

< |log(qi)-log(q^)| >l~Cy , C = — ~ 16.0 a s . (8) 

7T 

The transverse momenta are not ordered like in the Altarelli-Parisi evolution, but they perform a 
random walk. This means that there can be contributions to the ladder from configurations both 
with low and high transverse momenta. For the relatively small starting q^ we consider below, the 
momenta diffuse quite rapidly to larger values of q^_. 

• Rapidity distribution: 

Eq.(|^) can also be used to obtain the gluon density in the ladder at a given rapidity difference y. 
In Refs. |s|,|6| the inclusive cross section for the production of a gluon (gluon jet) with transverse 
momentum and rapidity Ya + y in the exchange of a BFKL Pomeron was derived for the case 
Y — y ^> 1 and i/> 1. Under these conditions the emitted gluon is not too close to the ends of the 
ladder. This cross section corresponds to fixing the momentum of a rung in the gluon ladder of Fig.l 
and can be described by Eq.(|J) for the parts of the ladder above and below the gluon in question. 
Upon integration over kj_, one can obtain the rapidity distribution p(y), the gluon density per unit 
rapidity interval of the emitted gluons at a given rapidity difference y, 

M -/ rfk ^^^^^ki)] 

- 1.6 a^/^pl . (9) 
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In the center of the spectrum one has p(y = Y/2) ~ 0.8 • aj 2 ■ VV. The BFKL results are valid 
only in multiregge kinematics, where < Sy > y ~ p^ 1 1. To fulfill this condition one should have 
a s Y <C 1.2a~ 2 . At the edge of the inclusive spectrum (y <C Y) one finds 

p(y)~ 1.6 a 3 J 2 Vv + D(a s ) , (10) 

where an unknown y independent function D which determines the density at small y, has been 
added. The rapidity interval between neighboring gluons < Sy > y — 0.6/(a« .^/y) becomes thus 
smaller with y towards the middle of the ladder. Note that the analytical result in Eq.(|Io|) can be 
simply understood. As the only dependence of the amplitude T n on Sy is due to the propagator, 
Eq.(S), this p(y) can be estimated as p(y) ~ 1/ < Sy > y ^ 2e(< >) ~ c?J 2 ^/y, where we used 
Eq.(£) to approximate < logq^ >~ y/a s y. 



III. MONTE CARLO APPROACH FOR THE POMERON 



To get a more detailed understanding of the hard Pomeron, a direct evaluation of the multidimensional 
integrals in Eq. (jl|) would be necessary. The expected gluon multiplicities for large total rapidity intervals 
are large, n ~ 10 4- 100, and thus the integrals of such high dimensionality that straightforward evaluation 
is impossible. We therefore have developed a Monte Carlo approach which enables us to deal explicitly with 
the multidimensional integrals. By keeping track of the kinematics, we can then examine to what extent 
the BFKL Pomeron indeed is made up from hard gluons and, for example, if the assumptions about the 
large rapidity intervals, Syi ^> 1, are justified throughout the ladder. An advantage of working numerically 
is that we can investigate what effect the use of a running coupling constant will have on the structure of 
the gluon ladder. We use for the running coupling constant 



1 

6 log(aTqi/A 2 ) 



'(<&) = ! TZ7T„ , „ 2 /A2 ^ . (11) 



i.e. the standard expression with A = 250 MeV, but extended by an additional parameter a that fixes the 
coupling constant at |qj_| = 0, thus regularizing the infrared singularity. The parameter a can be interpreted 
as the interaction with non-perturbative vacuum fluctuations which freezes the coupling constant at low 
momentum scales The value of a should be such that a s (0) ~ 0.4 -4 0.7. In our calculations we took 

a = 10.3, corresponding to a s (0) = 0.6 . 

As in the discussion in the last section, we will focus on the gluon subprocess indicated by the box 
in Fig.l. Therefore, we will generate for our Monte Carlo evaluation configurations {n,yi,]s.ij_} for the 
emitted gluons or, equivalently, for the intermediate gluons along the ladder, separated in rapidity by at 
least Sy m i n > 0, which is an input parameter. This is done in two steps. First, we generate ladders with a 
certain distribution in the momenta and in the number of rungs by an approximation method as described 
below. Each of these configurations is then used as the starting point for a Metropolis procedure to generate 
new configurations. This serves to correct for errors inherent in the approximation in the first step. This 
procedure leaves the distribution in n unchanged. 

In generating these multidimensional configurations, we make use of the specific dependence of the 
integrand on the kinematical variables, 

n 

\T n \ 2 ~Y[f(Sy i+1 ,c Li± ,c li+1± ) , (12) 

i=l 

3a, (qi , ) e^+i^q'+ii) 
f{Sy i+1 ,q A± ,q i+1± ) = = j- . (13) 
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Note that the external kinematics requires qu_ = and q n +i_L = qs_i_- Further we take for the gluon 
trajectory the usual form 

J (2tt) z (kj_ + /U )([k± - q±f + m ) 

but now with a running coupling constant and which at high q^_ for fixed coupling constant leads to the 
asymptotic form of Eq. (||) . 

The fact that Eq.@ is a product of functions of pairs of kinematical variables referring to subsequent 
rungs in the ladder allows one to generate the ensemble of configurations in a sequential fashion, for which 
we chose the von Neumann rejection method. It is essential that the factors f{5yi+\, qi_i_, qi+i_i_) are positive 
definite and thus can be interpreted as probabilities for the production of the (i + l)th rung with rapidity 
difference fo/i+i and linked to the former rung by a gluon with momentum qi + u_, when the previous 
rungs have already been generated. In fact, the conditional probability to produce the (i + l)th rung with 
(ch/i+i, qj+ix) as part of the whole ladder is given by 

f(Syi +1 ,qi±,q i+1 ±)-ip(Y-Y i+1 ,q i+1 ±) , (15) 
where Yi + \ = Yl % k=i is the rapidity interval between the 'zeroth' {i.e. particle A) and (i + l)th rung and 



oo . n+l / n+1 \ 

i>(Y-Y i+ i,q i+1 x)= X) / II dyi d2 <ii±f( s yi, c ii-ii-,<ii±-) s {y-Yi+i - 5 yi) 

n=i+l l=i+2 \ l=i+2 J 



(16) 



is the total probability of configurations 'after' this (i + l)th rung. Since the integrand in if) is given by a 
similar factorized form as Eq. (|l2|) , we can use Eq.(H|) to see that at very large Y — Y i+ i 

^(y-F i+1 ,q i+lx )~ex P ^- m y_ Yi+i) j , (17) 

where qsx is the transverse momentum of the final gluon in the ladder which is taken to be fixed and not 
too large. So for Y — Yi+% ^> Yi+\ the function ip is constant for values of logq^ +1 ~ yi^+i which 
we expect to be essential in the generation of the quantities (dyi+i, q^+i) at the (i + l)th rung. Therefore 
when generating the initial part of the ladder, it is a good approximation to use only the first factor 
/(<%i+ij qiJ.) qi+iJ_) in Eq. (|i~5|) . Starting with the first rung, for which we take |qA_i_| = -R^ 1 ' we can 
therefore generate subsequently for every rung the variables (Syi,<lii_). We repeat this until Y n+ \ ~ Y , 
which in our approach then fixes the gluon multiplicity n in this ladder. We stress that we will only use the 
first part of the ladder for the results below. 

It should be noted that when using the weight functions generated with the von Neumann rejection 
method, we have to introduce an ultraviolet cut in k^. This is due to neglecting ip in Eq.(|l5|). The 
dependence on this cutoff parameter Q m ax is of type log log ((^^//i 2 ) and thus weak. Since we fix the end 
momentum of the ladder, |qsjj ~ -R^ 1 , the dependence on Q m ax vanishes. 

For the Metropolis routine, we use the configurations obtained by the von Neumann method as 
starting configurations to improve the sampling in the phase space, producing a series of new configurations 
{n, y[, ]s.'i±}. Note that the number of rungs n is not changed. A more detailed discussion will be published 
elsewhere. 

To calculate transverse momenta distributions, mean rapidity intervals etc., we take a very large total 
rapidity interval Y and consider only a part of the ladder for which y <C Y in accordance to our approx- 
imation -0 ~ constant. That this yields predictions independent of Y and other quantities characterizing 
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the particle B at the end of the gluon chain can also be seen, for example, from the asymptotic expression 
for the inclusive cross section, Eq.(^|). So one should realize that using our method we can investigate only 
distributions in the growing part of the gluon spectrum. 

For the dependence of the total cross section on the energy s we have to proceed differently. In order 
to obtain the Pomeron intercept it is sufficient for the total cross section to integrate ( |l2| ) over the full 
phase space with logarithmic precision. This was done in two ways. In the first method, we used the ratio 
of accepted to total events in the von Neumann procedure at different values of Y ~ logs. Fitting this 
to an exponential function, we obtain the exponent A in Eq.(Q). In the second method, we first calculate 
for a given energy s (or rapidity interval Y) the mean multiplicity of the ladder and then approximate 
the Pomeron by a ladder containing that many rungs. In the corresponding amplitude, Eq.(|12|), we then 
replace the transverse momenta of the previous rungs in the products of the functions f(5yi+i, qi_i_, q^+i^) 
by the already determined mean value for that rung. This approximation enables us to compute all integrals 
needed for the total cross section separately. Performing this again for different values of Y the exponent 
A can be extracted. Both methods were found to agree reasonably well. 



IV. RESULTS AND DISCUSSION 



We first show what our Monte Carlo calculations yield for the properties of the gluon ladder discussed 
in Section 2 and compare them to the predictions of the BFKL theory in the asymptotic limit. In all our 
calculations presented in this section we chose for the "constituent" gluon mass /i one half of the gluonium 
mass: /i ~ 0.5 GeV. We have checked that changes by a factor of two in this parameter had a negligible 
effect on the results below. Furthermore, we established in all cases that there is indeed no dependence on 
the ultraviolet cutoff introduced in the von Neumann step. 

First we examined the exponent A, which fixes the Pomeron intercept. We performed calculations 
with a fixed coupling constant in the range 0.01-0.20 and for a total rapidity interval Y = 50. For the 
lowest values of a s this yields rather short ladders (n ~ 6) for which our method of generating rungs in the 
ladder is less reliable. For a s = 0.05, we deal with ladders of total length n ~ 20 and obtain A = 0.11, close 
to the BFKL result of Eq.(Q) which gives 0.13. Taking a s — 0.1, we get numerically 0.27, compared to the 
BFKL value of 0.27. This reasonable agreement indicates that our numerical approach is quite reliable. 

To use our samples to extract the gluon transverse momentum and rapidity distributions, we should 
restrict ourselves to the first part of the ladder, i.e. up to rapidity intervals of y ~ 20 in our example. For 
the final gluon we took |qs_i_| = Rg 1 ~ . In this special case with similar transverse momenta for both 
ends the last part of the ladder is identical to the first part and, aside from some uncertainty in the middle of 
the ladder, the complete ladderstructure can be obtained. Results for the mean logarithm of the transverse 
momenta for an initial gluon momentum of Iq^iJ = R^ 1 = 10 fm _1 and for two values of the coupling 
constant a s are shown in Fig. 2. Indeed we see that < | log(q^/q^ ± )| > 2 follows the behavior predicted 
by the BFKL theory, Eq.(||): a linear rise with the rapidity difference y. The slopes also agree reasonably 
with the BFKL prediction. For a s = 0.05 and a s — 0.1 we find respectively 0.31 and 0.51 compared to 
the analytic BFKL values of 0.8 and 1.6. The gluon density p(y) is shown in Fig. 3. Here our results also 
reproduce the BFKL predictions well, both the y dependence and the coefficients; we fitted both curves to 
the function p(y) = Cq + C\^fy and found for C\ a value of 2.2 10 -2 for a s = 0.05 and 5.8 10~ 2 for a s = 0.1 
wheras BFKL in the asymptotic limit predicts 1.8 10~ 2 and 5.1 10~ 2 , respectively. 

This good agreement with the main BFKL features confirms the validity of our method. We now 
use it to examine several aspects that go beyond the original framework of the BFKL Pomeron. The first 
concerns the effect of a running coupling constant. We repeated the calculation for the total cross section 
with the a s given in Eq. ( ji"l| ) • If we try to fit the energy dependence according to Eq. (Q) , the effective value 
for A clearly decreases compared to that obtained with a fixed coupling constant a s = a s (q\j_), but our 
method of calculating is not reliable enough to give its precise value . This decrease can be understood 
from Fig. 2, where we also show a curve for a running coupling constant. We observe that q^ continues 
to grow as we go down the ladder (just as for fixed a s ); this is due to our choice of q 2 A± which leads to 
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a rapid diffusion from log(q^ ± ) mainly to larger values of log(q^). This leads to contributions with lower 
coupling constants, resulting in a lower A. The effect of a running coupling constant on the gluon density 
is shown in Fig. 3. While there are some changes, it has no essential influence on the qualitative features of 
the distribution of gluons in the hard Pomeron. There have been previous studies of the effect of a running 
coupling constant f|,||,D ■ They showed the decrease of A and its influence on the character of the Pomeron 
singularity. 

Another important aspect we can investigate is the rapidity interval, 5y, between successive rungs of 
the gluon ladder. What we show in Fig. 4 is the average rapidity interval as function of the rung number in 
the ladder. Our results show a rapid drop of < Sy > for the first rungs, going over in a much slower decrease 
which continues to the end of the ladder. In our calculations we imposed a lower bound of Sy rn i n — 0, since 
the leading contribution is assumed to come from high 5y configurations. Fig. 4 makes clear how important 
the value of this cut is: a major part of the ladder consists of rungs with only a relatively small rapidity 
separation, typically of order 1.5 for a s = 0.1. Thus had we chosen a large 5y m i n , entirely different results 
would have been obtained. The use of a running coupling constant leads to the same observations. For a 
smaller fixed coupling constant {e.g. a s ~ 0.05) we find rapidity intervals mainly around 3.0, which means 
imposing 5y m i n > would have a smaller effect on the ladder structure. From this analysis we see that 
configurations with small rapidity intervals in the ladder play an important role and can influence the results 
for larger values of the fixed coupling constant and for the running one. 

Our Monte Carlo study of the hard Pomeron in terms of explicit calculations of gluon ladders repro- 
duced the main predictions of the BFKL Pomeron over a wider range of kinematical conditions. This seems 
to confirm the validity of the kinematical assumptions that go into the BFKL derivation. However, when 
one looks in more detail at the gluon ladder, which can be done with our numerical method, it becomes clear 
that the multiregge kinematics are not applicable anymore towards the middle of the ladder, especially for 
not too small values of a s . In a forthcoming publication, we will therefore study the validity of the BFKL 
Pomeron in more detail for different kinematical conditions and the importance of higher order corrections 
through numerical simulations. Theoretical progress has been made in computing the next-to-leading order 
terms [ ^0[ , but the expressions are rather complicated and there is presently little hope that higher order 
terms can be computed exactly. More intuitive pictures for the hard Pomeron such as our direct ladder 
calculations or the color dipole model p"l| ] can show the way towards a better understanding of the nature 
of the "hard" Pomeron. 
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FIGURE CAPTIONS 



Fig.l The contribution of a gluonic ladder with n rungs to the Pomeron exchange. The Lipatov vertices are denoted 
by the dots and the reggeized gluon propagators by the thick, vertical gluon lines. 

Fig. 2 The square of the mean for the absolute value of the logarithm of the gluon transverse momentum along 
the ladder , | log(q5_/qA±)|j as function of rapidity interval y for different coupling constants; the dashed line for 
a s = 0.05, the solid line for a s — 0.1 and the dot-dashed line for a running coupling constant. 

Fig. 3 The gluon density, p, as a function of rapidity interval y for different coupling constants; the dashed line for 
a s = 0.05, the solid line for a s = 0.1 and the dot-dashed line for a running coupling constant. 

Fig. 4 The average rapidity interval between two neigboring rungs as function of the rung number N for different 
coupling constants; the dashed line for a s — 0.05, the solid line for a s —0.1 and the dot-dashed line for a running 
coupling constant. 
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